arXiv:1502.02752vl [nlin.SI] 10 Feb 2015 


On an integrable system related to the relativistic 
Toda lattice -Backlund transformation and integrable 

discretization 

Luc VinetJ, Guo-Fu Yu* *and Ying-Nan Zhang^ 

* Centre de Recherches Mathematiques, Universite de Montreal, 

C.P.6128, Centre-ville Station, Montreal, Quebec, H3C 3J7, Canada 
•^Department of Mathematics, Shanghai Jiao Tong University, 

Shanghai 200240, P.R. China 

§LSEC, Institute of Computational Mathematics and Scientific Engineering Computing, 
Academy of Mathematics and System Sciences, 

Chinese Academy of Sciences, P.O. Box 2719, Beijing 100080, P.R. China 
•^Graduate School of the Chinese Academy of Sciences, Beijing, P.R. China 


Abstract 

We study an integrable system related to the relativistic Toda lattice. The bilinear representation of 
this lattice is given and the Backulund transformation obtained. A fully discrete version is also introduced 
with its bilinear Backlund transformation and Lax pair. One-soliton solution of the discrete system is 
presented by use of Backlund transformation. 
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1 Introduction 

In terms of the canonical variables {q n ,Qn}n£Z, the Hamiltonian for the infinite relativistic Toda lattice 
(RTL) is given by M 

H{q, 0) = ^{exp(0„)[l + exp(g n _i - q n )} 1/2 x [1 + exp(g n - q n+1 )} 1/2 - 2}. (1.1) 

n(zZ 


The equations of motion are hence 

d _ dH n _ 
dt qn ~ 8Q n ~ n ’ 

d _ dH n _ 1 . . 1 

,, F'n. — ^ — n ^n—l(^n ®n— l) r% ^ , n\Pn "F ^n+1/5 

at oq n 2 2 

with 

K = exp(0„)[l +expQ„_i - q n )] 1/2 x [1 +exp(g n - q n +i)], 
_ exp (q n - q n+ i) 

1 + exp(g n - q n +i)' 

* Corresponding author: gfyu@sjtu.edu.cn 


( 1 . 2 ) 

( 1 . 3 ) 


( 1 . 4 ) 

( 1 . 5 ) 
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The RT equation 


.. _/, , 1 • V1 . 1 • x exp(g„_i - q n ) 

Qn —(1 + -9n-l)(l + -Qn) . n / 2 \ - 1 -\ 

c c 1 + (l/c 2 )exp(g n _i - q n ) 

n . 1 • V i , 1 • \ exp(g„ - q n +i) 

- (1 + -g n )(l + —Qn+l) 1 , / 2 s- 7 - T, 

c c 1 + (l/c 2 )exp(g„ - q n+1 ) 


( 1 . 6 ) 


where q n is the coordinates of n —th lattice point, and means the differentiation with respect to time t and 
c is the light speed, was introduced and studied by Ruijsenaars Q. The evolution equations of {a n ,b n } 
assume the form 


7, — ^n(^n—l^n—1 ®n^n+ 1)5 

at 

(1.7) 

= ®n(l ^n){pn ^n+l)* 

(1.8) 

By a further transformation 



(1.9) 

the equations (jl.7p~(|1.8|) are rewritten as 


"TT^n = UniVn— 1 ^n)i 

at 

(1.10) 

d . v 

"TT^n — 'Vny'Vn— 1 ^n+l “1” ^n+lj- 

at 

(1.11) 


We refer to dl.10ll - dl.lll) as the RTL. HI id- 

The integrable lattice equations are related to the following discrete spectrum problem and time evolution 
equation 


Eip n = Un(u, X)ip. 
<Pn,t = V n (u,X)(p. 


n-i 

n? 


( 1 . 12 ) 

(1.13) 


where u is the potential function, A the spectral parameter and E the shift operator defined by Ef n = f n +i ■ 
The integrability condition between (11.1211 and (11.131) leads to the integrable lattice system 


U n ,t + U n V n — V n+1 U n — 0 , 


(1.14) 


when we take u = (u n , v n ) T , and 


U n = 

V n = 


f \ 2 +u n A\ 
\ Xv n 0 ) 


A 2 /2 - v n -i 
Xv n —i 


— A 2 /2 -Un - 1 


(1.15) 

(1.16) 


The compatibility condition (11.141) entails the RTL equation dl.lOt - dl.llI) . So (11.121) and (11.131) with (11.151) - 
(11.161) constitute a Lax pair for the RTL (ll.10D - dl.llD . Other Lax presentation for dl.lOD - dl.llI) can be found 
inp-Q. The bilinear form and Casorati determinant solution for the RTL (EH) were given in Q. It is 
also worth pointing out that the RTL dl.lO|) ~ dl.lll) can be related to the Laurent bi-orthogonal polynomials 
HQ The discretization of the RTL dl.6D was first performed in Hi using the Hamiltonian method. The 
discrete-time relativistic Toda lattice (dRTL) equation was first proposed by Suris 


A exp(g^ +1 ~ q^) - 1 = 1 + g 2 exp(q t n _ 1 -q t n )l + ( g 2 /S) cxp(g* - g*+\) 
8 exp (qi - gn _1 ) - 1 1 + 9 2 exp(g*. - g‘ +1 ) 1 + ( g 2 /5 ) exp(g£l\ - g£) 


The Lagrangian form of dRTL was presented in [18| . The Casorati determinant solution for the discrete RT 
lattice (II. 171) was given in 141 and the elliptic solutions in Hi- 
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Although the discrete version of the RTL (11.61) is already known, as far as we know, the discrete analogue 
of the RTL related system 11.101) - (11. 1 ID has not been given. The purpose of this paper is to propose a direct 
discrete version of the RTL (I1.10l) - (ll.lll) using Hirota bilinear method. It is nontrivial and of considerable 
interest to find integrable discretizations for integrable equations. Attention is being paid to the problem of 
integrable discretizations of integrablc systems. (Seee.g., [16h2C| and references therein). Various approaches 
to the problem of integrable discretization are currently available. One of them is Hirota’s bilinear method 
[Umi, which is based on gauge invariance and soliton solutions. Here, we focus on a discretization process 
such that the resulting discrete bilinear equations have Backlund transformations (BTs). As a bonus, we 
can usually derive Lax pairs for the resulting discrete equations. This method has been successfully applied 
to the discretization of (2 + l)-dimensional sinh-Gordon equation 27f and of the two dimensional Leznov 
lattice equation 28j. Based on bilinear forms and determinant structure of solutions, Hirota’s discretization 
method has also been developed to construct discrete versions of the Camassa-Holm equation 129] and the 
short pulse equation |30]. 

The content of the paper is organized as follows. In section 2, we give a bilinear form of the RTL. In 
section 3, a BT and Lax pair are exhibited. In section 4 an integrable discretization of the RTL is presented 
and the integrability is made manifest by the corresponding BT. Section 5 is devoted to conclusion and 
discussions. 


2 Bilinear form for the RTL 


Performing the change of the dependent variable 

d , F n 

v ’‘ = dt l "G7 l + h 


d , G n _i 

n = dt ~F\T + *’ 


the system (11.1011 - (11.1111 is transformed into the following bilinear form 

DtGn— 1 • F n + F n G n — \ — C\F n — 1 G n , 
DtF n • Gn F n Gn — C2F n —\Gn+1 3 

which can also be written as follows in terms of the Hirota bilinear differential: 


(D t e 
(D t + 1 - c 2 e 


— e 2 "" + cie 2 
-d„ 


')F,G = 0, 


')F.G = 0. 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 


(2.4) 

(2.5) 


Here ci,c 2 are arbitrary constants. In the following we shall take c\ = c 2 = 1 for simplicity. The Hirota 
bilinear differential operator D^ and the bilinear difference operator exp (SD n ) are respectively defined [31] 

by 


D^a»b= 


a{t)b(t')\ t '=u 


d__d_ 

^ dt dt' 

exp(SD n )a(n) • b(n) = a{n + 6)b(n — 6). 

It is well known that nonlinear integrable equations share many common features, among which, the BT 


and their associated nonlinear superposition formulae [32 


3 Bilinear Backlund transformation for the RTL 

For the sake of convenience, we introduce an additional discrete variable m and set 

F = f m+ 1, G = / m _i 

We shall denote F[t , n) by F n and f(t, n, m) by /„ im for simplicity. With the help of (13.11) . the eqs. 
reduce to 


(3.1) 

3, m 


(D t e- 
(■ D t e - 


+ D m Dn + Dm D m -D n x 

2 -e 2 +e 2 )/./ = 0, 


2 — 2 


■--D r , 


)/•/ = o. 
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(3.2) 

(3.3) 























Proposition 3.1 The bilinear equations t3.2\) . \3.S\) have the Backlund transformation 

(D t -Xe~ D " +fi)f.g = 0, (3.4) 

, Dn-Dm Dm~P>n D n +D m . 

(e 2 — e 2 — Ae 2 )f*g = 0, (3.5) 

(. Xe~ Dn ~~^ +e~~^ + r je~^ L )f»g = 0, (3.6) 


where A, /r and 7 are arbitrary constants. 

Proof. Let / be a solution of equations (I3.2M3.3II . If it can be shown that g given by (13.4M3.6I1 satisfies 

Dn + Dm. Dn + Dm Dm-Dn 

Pi = (D t e 2 -e 2 +e 2 )g»g = 0, 

P 2 = (D t e^P + e^~ - e~^ L - Dn )g»g = 0, 

then the equations (13.4M3.6I) define a BT for (13.21) - (13.31) . 

In fact, using (13.41) - (13.61) and the bilinear identities, it can be seen that 


-[e— f.f]P 2 

= [( D t e^ +e^P - e^P~ Dn )f •/] [e^P g»g] - [(D t e^ + e^~ - e^~~ Dn )g . g}[e^ f • f] 
= 2sinh {^)[D t f ,g\.(fg) - 2 sinh( ^) (e ° n 2 ° m / . g) . (e D "* 2 / . g) 

= 2 sinh( ^ L )[D t f.g\.(fg) + 2A sinh( ^-){fg)^(e~ Dn f.g) 

= 2sinh(^)[(A - Xe~ Dn )f »g\ • (fg) 

= 0 


r Dm+Dn 

— [ e 5 f • f]Pi 

r _ Dn+Drr, Dn+Dm D m ~Dn , n D m + D n 

= [(Ae 2 -e 2 +e 2 )/./][e 2 g»g] 


Dn+Dm Dn+Dn 

- [(Ae 2 - e 2 


- e 

, d 


)ff«5][e Dm ^"/«/] 


= 2 sinh( Dm + D n ) (A/ . g)»(fg)- 2 sinh( ) (e 2 ‘f.g).(e 2 ' f»g) 

= 2 sinh( Dm + ^ )(Ae~ r> ’ 1 / « g). (/g) - 2sinh(^)(e ; ^ IL /.g). (e’^/.g) 


= 2sinh(^)(e^/.g).[(-Ae 
= 0. 


-D n ~- 


2 — e 2 


')/» 3 ] 


This completes the proof of proposition 2.1. 

From Proposition 2.1 it can be deduced that the bilinear RTL (I3.2M3.3I) has the following Backlund 
transformation: 


D t F.F' - AA-i^+i + /uFF' = 0, 

(3.7) 

AG.G' - AGn-iG^i + nGG' = 0, 

(3.8) 

G„ +1 A - FG; +1 - AGF' +1 = 0, 

(3.9) 

AG n _iF^ +1 + GF' + jFG r = 0, 

(3.10) 

F'= g m+ 1 , G' = g m _ 1 . 

(3.11) 


where F' and G' are defined by 

Starting from the bilinear BT (13.711 - (13.101) . we can derive a Lax pair for the system (11.101) - (11 .111) . First, set 


a, _ , T , _ Gn p _. . F' 

^ra — p, 1 ^ ! Pn — In p/ , Q n ~ l n 


G' 


F' 


G' 


4 




























in (13.7I) - (13.10I) . Eqs (13.7D - (13.101) are then transformed into 


V - Ae Pn+ ^ n_Pn+1_< ^ n+1 4' n _ 1 + ^ „ = 0, 

(3.12) 

■n,t - \ e P "-+Qn-l- p n + l-QnQ n _ 1 + ^ = Q; 

(3.13) 

'n+1 ~ ~ Xe Qn+1 ~ Qn ^ n = 0, 

(3.14) 

'n + 7$n + Ae Pn_P " +1 4'„_i = 0. 

(3.15) 


From the compatibility of above equations we now obtain 


dt 

d 


rP n - e Q ^- Q - = —P n+1 - e Q -- Q 


dt 


Q n _ gPn-Pn + l+Qn-l-Qn + l _ __ g-Pn+1 ~Pn+2+Qn ~Qn+2 , 

dz n di n 


so we have 


— Pn = e 
dt 


— e Qn-l—Qn 


+ Cl, 


T .Qn = e P "- P "+ 1+C? "- 1 - Q '*+ 1 + 02 , 

at 


(3.16) 

(3.17) 

(3.18) 

(3.19) 


Here ci, C2 are two arbitrary constants. We now introduce the new potentials u n = P n .t — C\,v n = Q n ,t — C2 
and differentiate Eqs. (13.181) and (13.191) to find 


dt 

d_ 

dt 


u n = e Qn ~ 1 Qn {v n -i - v„) = u n (v n -1 - Un), 


_ pf’n-f’n+l+Qn-l-Q? 


Vn = e 


+ 1 (^n %+l "I" ^n-1 ^n+l) 


— ^n(^n ^n+1 H - ^n—1 ^n+l)j 


(3.20) 


(3.21) 


which is the original nonlinear RT lattice (11.101) and (11.111) . So (I3.12D - (I3.15D provide a Lax pair for RTL 

ehd-eid- 

Using (13.71) - (13. 10D . we can obtain a nontrivial solution from the trivial solution. It is easy to check that 
F' = G' = 1 is a solution of (12 .21) - (12.31) . Substituting this trivial solution into the BT (13.71) - (13. 101) . we have 


Ft — \F n _i + /. iF n — 0, (3.22) 

Gt — AGn-i + IJ-G n = 0, (3.23) 

G n+ i -F n - AG„ = 0, (3.24) 

AGn-i + G n + ”fF n = 0, (3.25) 

When we take A = 2, /z = 3, 7 = 3 we can obtain the following solution of m-m 

F n = (2/3)" +1 — 2(2/3)" — e -t , G„ = (2/3) n + e _t . (3.26) 


In Fig. 1, we plot the one-soliton solution for u and v respectively. 


4 Integrable full-discrete version of the RTL 

We propose the following bilinear equations 

r 1 , . Dm + Dn D m + £) n , ^ Dm — D n , 

[_ sinh(eD t )e- 5 — -e- 5 — +eD t+& — 5 - ejDt ]/-/ = 0, (4.1) 

[—isinh (eD t )e~~P + e^~^~ eDt — e~? L ~ Dn ~ eDt ]f • f = 0. (4.2) 

In the continuum limit as e —> 0, the system (ED and (14.21) reduces to the lattice (13.21) and (13.31) . The system 
ED and (14.21) therefore provides a fully discrete version of (13.21) and (13.31) . 
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Fig. 1: (Color online) The one-soliton solution for (11. 101) - (11.Ill) : n = 4. 


Note that 

F = f m+ G = f m _x, 

and denote the variable t as fc, Eqs. (IXTl)-(lTa can be rewritten as 

/ ^ i \/^ik+l j7>k—l , ^ ryk—1 ZT'fc+l /^fc+1 zt'A;— 1 n 

V 2 e ^n+1 ' 9 6 ^ n ^+1 ' ^n+l-^n — 

(1 + lJGjT 1 ^ 1 - - G k n %\F*z{ = 0. 


With the dependent variable transformation 


s~ik 

ui = n ~ l 


rpk 

V k = — 
pk ’ n ’ 

-t M n 


the bilinear equations (14.4D - (14.5I) are put in the following nonlinear form: 

( _ \ _T)w fe + 1 ,, fc + 1 -4- ii k ~ 1 v k ~ 1 4 _ ^-ii k ~ 1 v k+1 — 0 

(— 4 - lb; fc + 1 7 / fe+ 1 i ; fc+1 _ L v k-KMl,.k+l _ k-1 k-1 k-1 _ q 

\r, e + ij« n “n+lCi+l t) e U n “ra+lTn+l U n u n-l u n — u - 

Set t = ke. By using the Taylor formula, 

u k n +1 = u n {{k + l)e) = u n {t + e) = u n {t) + u' n (t)e + 0{e 2 ) 

and after careful calculations, we finally obtain that as e —> 0 the continuum limits of the fully 
(1377b and (14.811 are given by 

, / Wn-l , 

u„ = u n |-1 

v n 

'Un'Vn—1 
v ^n+l^n +1 

Setting d n = (lnu n y + 1 and c n = (lnu„)' + 1, then we get 

d n — dniCn—l O n ), 

C n — Cn(Cn — l Cn4~ 1 T ^n+l); 


-1 . 


(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 
discrete Eqs. 

(4.10) 

(4.11) 


(4.12) 

(4.13) 
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which is just the RT system (11.101) and (11.111) . Therefore the system (14.71) and (14.81) provides a fully discrete 
version of the potential RT system. 

When we take the dependent variable transformation 


Tpk— 1/^fc+l 
Tjk _ 1 n ^n— 1 
Uqn ~ k+l n k-l 


Ff +1 G: 


we can deduce from (mu-el that 

1 


(1 

(1 


< +1 

Un~ l 1 


v k = 


pk+lr<k-l fc-pi 

J- n n '-'n 


rpk — 1 s^ik -\-1 

F n 


k -2 


Vn 1 ’ 


— ){U k+1 - U k ~ l ) = n ~ 1 

2e >{ n n ’ v k + 2 


1 


u n—l 


+ 7 r.Wn -V* ) = 


2 e 


One can check that in the continuum limit 


# fc +2 fc +2 


k-2 k-2 
l n u n— 1 


*n+l u n+l a n+l y n+l 

0, Eqs. (14.15l) - (|4.16l) reduce to the RTL equation 


(4.14) 


(4.15) 

(4.16) 


U' n = U n {V n - X -V n ) t (4.17) 

Vn = V n {V n _i - K +1 + u n - U n+1 ). (4.18) 


Hence equations (14.14l) - (14.16l) give a fully discrete version of the RT system (11.101) and (11.111) . 

Remark 4.1 We remark here that the full discrete analogue of RTL obtained here is different from (11.171) . 
Suris obtained dRTL by using Hamiltonian method, while here we start from Hirota bilinear method. What 
interesting for integrable discretization is that one integrable system may have some different integrable 
discrete analogues. Among these integrable discrete versions, some are linked each other by appropriate 
transformations. Others may have no direct relation. In our case, we have not found the direct link between 
our full discrete RTL to (11.171) . 

For simplicity we take e = 1 in the following. In this case, the bilinear system m-m and nonlinear 
equations (14.7I) - (14.8[) become 


and 


‘--D k 


e 2 


-D n —Dk 


3 

~~ 2 6 
1 

+ 2 6 " 


L 


1 Dm+Dn ^ , 

776 -5- ° k )f • f = 0, 


2 6 


~^~ Dk )f-f = 0, 


(4.19) 

(4.20) 


ut^vt- 1 - \u k n \\v k n X 1 + = 0 , 

1 

2 l 


' t, n+1 u n 
.k 




*n -\-1 u n+l 


1 v k ~ 1 v k ~ 1 4- A v k -i ll k + 1 v k + 1 _ -v k+1 ii k+1 v k+1 - 0 

""-V" I ~ U n. a n+l u n+l u n u n+1 u n+1 —v. 


(4.21) 

(4.22) 


respectively. 

We will show that (14.21l) - (14.22l) is integrable by exhibiting its BT and Lax pair. Concerning the bilinear 
equations (14.191) - (14.201) . we have the following result: 

Proposition 4.1 The bilinear equations \4.19ty~lf.2Cf) have the Backlund transformation 


(e + ye - 2/je”” D ")f.g = 0, 
(e^ L+Dk + \e~^ L ^ Dk + ge~^ L+Dk )f »g = 0, 

Dm.+Pn. 2 P m -p„. Dm+Dn . 

(e 2 + -pe 2 + ae 2 )f.g = 0 . 


(4.23) 

(4.24) 

(4.25) 


where A, fi, 7 and a are arbitrary constants. 
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Proof. Let / be a solution of equations (14.19l) - (14.20l) . If it can be shown that g given by (I4.23I) - (I4.25I) 
satisfies 


Pi = (e 
P‘2 = (< 


^-D n -D k 


3 Dm+D 

~ 2 6 
1 


-+D k 


1 D m +D 
2 6 


+ 


-i^+D k _ ^ 
2 


Dk )g*g = o, 
Dk )g*g = o, 


then the equations (14.231) - (14.251) provide a BT for (I4.19I) - (I4.20I) . 

In fact, by using fl4.23D - fl4.25D and the bilinear identities, one can see that 


'--D k 


f»f]Pi 

Dm—Dfi 3 Dm~\~D n . t—\ Dm, + Dn ir Dm, + Dn t-\ 

— [(e 5 Dk - 2 +Dk + -e 2 fc )/-/][e * Dk g»g} 


- - -e 


3 Dm + Dr, 


-\-D k 


1 Dm + Dn 


~D k 


)g»g][e 


Dm+Dn 


~D k 


/•/] 


2 2 6 

= 2sinh (^Y-)(e~^ L ~ Dk f . g). {e^P +Dk f .g) - 3 sinh(D fc )(e _ / . 3 ) . (e Dm t Dn f . g ) 

= -2g sinh(^) (e~ ^ ~ Dk f . g) . (e~ ^ +Dk / . g ) - 3 sinh [D k ){e~ Dn \° m f.g). / . g) 

= 2g sinh(D fe ) [e D ” 2 f • g] . (e~ Dm t Drl f . g ) - 3sinh(D fc )(e" Dn \° m f. g ) . ( e Dm i° n f, g ) 

Dm + Dn D m + D n 2 D m -D n 

=-3sinh(D fc )(e 2 f.g).[{e 2 + 2 )f.g\ 

= 0 


and similarly that 




= l(e^- D "~ Dk + ±e~^ +Dk - “e-^~ Dk 


1 Dr, 

~2 


3 D n 

2 ( 


)/•/][« 


e 2 


■-D k 


g»g\ 


-[(« 


" Dn D ' c +^ e ° m+Dfc -^e Dfc )5»<?][e' 


*/•/] 


= 2sinh( ' Pm 2 Dn )(e Dk f»g) • (e^ 2 " +Dk f»g) + sinh(D fe )(e f • g) • f .g) 

= -2/rsinh( £>m ? Dn )(e~^ L ^ Dk f .g). (e~^ +Dk f . g) + sinh(H fc )(e~^/.g).(e^ L /.g) 

= 2/r sinh(D fc )(e^ _Dn / • g).{e~^ L f .g) + sinh(D fc ) (e - ^ - / .g).(e^f.g) 

= sinh (H fc )(e _ ^/.g).[(e^ i - 2ge^~ Dn )f .g] 

= 0 . 


This completes the proof of proposition 4.1. Again we can obtain from Proposition 3.1 that the bilinear fully 


discrete RT lattice fl4.4D - fl4.5D (e = 1) has the following BT: 

F n G' n + ^G n F' n - 2 M F n _ 1 G( l+1 = 0, (4.26) 

F n +iG' n + ^gF n G' n+ 1 + aG n F^ +1 = 0, (4.27) 

Fn+l,k+l P n ,k-1 F ^Pn,k— l-^n+ljfc+i F gPn,k+lF n+lk _i = 0, (4.28) 

G n +i,k+iG' n j._ 1 + \G n ^-iG' n+l k+1 -)- gG ny k+iG , n+l k _ 1 = 0. (4.29) 


Using (14.2611 - (14.291) . we can obtain a nontrivial solution from the trivial solution. It is easy to check that 



















































F' = G' = 1 is a solution of (14.4I) - (14.5I) . Substituting this trivial solution into the BT (I4.26I) - (I4.29I) . we have 


F n + 'yGn — 2^iF n -i — 0, 

2 

F n + 1 + —fj,F n + aG n = 0, 

F^+l+XF^ + ^+^O, 

G^W + XG^ + ^+^O. 

When we take a = II 7 ,fi= 1, A = —100, we get the following solution of (|4.30[) - (|4.33p 


22 

F„= Cl 5 fc 3" + c 2 (yi2) fe (-)", 


G n = 2 


22 , 


ci5 fe 3 "- 1 + c 2 (\/l2) fc (7) 


n— 1 


22 , 


-ci5 fe 3 n -c 2 (\/l2) fc (—)". 


(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 

(4.35) 


ci and c 2 are arbitrary constants. In Fig. 2, we plot the one-soliton solution for U k and V k by use of (14.1411 
with parameters selected as ci = 2, c 2 = 1 and k = 5. 




Fig. 2: (Color online) The one-soliton solution for full-discrete system (I4.14D - (I4.16I) . 


In the following, we derive a Lax pair for the fully discrete RTL (I4.21D - (I4.22D . Let 


F' 

n 

G' 


GL 


F’ 


F _ G T 

— = <F — = ip 

F' ’ G' 


From Eqs. (14.2611 - (14.2911 we get 


$ n u n+1 v n+1 + 7^„u n+1 v n +i ~ 2/J,$ n _ 1 U n V n _ 1 = 0, 


<& n+ 1 V n+1 + a^nVn+l + nVn = 0, 
$ k+ \u k ~\v k ~ 1 

^ n+l u n+l u n 




fc+1 fc-lfc-l 


X^ k ~ 1 u k ~ 1 v k ~ 1 - 

A ^n Uj n+l u n 
\ VTrfc — 1 „ , fc — 1 „, fc — 1 

- Aw„ u n+1 v n+1 


^n +1 U k +\v k+1 = 0 , 


' n+l u n+l u n+l T ''*n u n+l u n+l " r “n+1 u ra+l 

In order to obtain a Lax pair for (14.211) and (I4.22L we first use (14.3811 to express ’F in (14.401) in 
and then to write $ n+2 and 4> ra+ i as expressions of 4>„. After some calculations we find 


fc +1 fc +1 fc +1 


U 


V„ 


= 0 . 


-2v k ~ 1 v k ~ 1 v k ~ 1 v k+1 +- 2v k ~ 1 ii k ~ 1 v k ~ 1 v k+1 
^ u n+l u n+2 u n+l u n+l ^ zu n +l u n+2 u n u n+2 

_q 7 A+i 7 .k— l..fe+i i o 7 A—l„,fc+i„,fc+i„,fc+i 

ou 'n+l“n+2 ( V|-l Ci+2 ' ou » 


n+l u n+2 ^ra+l y n+2 






= 0 , 


(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 
terms of $ 

(4.41) 


or 


za n+2^n+l oa n+2 t; n+2 
fc— 1 fc +1 
u n+2 v n+2 


2u\ 


k-1 k-1 

n+l u n 


— 3 U[ 


k-\-l fc+1 
n+1 u n +1 


fc -1 fc +1 
^n+l^n +1 


2u k ~ 1 v k ~ 1 — 3?y fe+1 ?; fc+1 
= (E n - 1) W+1 , 7 n+1 " +1 = 0. 




fc— 1 fc +1 
n+1 u n -\-1 


(4.42) 
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It follows that 


or equivalently, 


2 n k ~ 1 v k ~ 1 — Hn k+1 v k+1 

zu ra+ l u n _ ou n+l u n+l 

fc -1 fc +1 

“n+l+ri+1 


2v k ~ 1 v k ~ 1 - 1v k+1 v k+1 

zu n+l u n ou n+l u n+l 


rl ,k— l,,fe+l — n 

cu n+l U n+l — U ' 


Setting u k 


c k / 2 u k ,v k —> c k / 2 v k , we have 


2v k ~ 1 v k ~ 1 — 4w fe+1 « fc+1 

zu n+l u n ou n+l u n+l 


fc - 1 ,, fc +1 _ 




= 0 , 


■‘'ra+l+ra+l 

which coincides with (14.211) . Similarly, combining (14.471) and (14.491) . we get 


- 2 u! 


fc —1 1 +fc fc —1 1 +fc fc —1 


n+l U n+l u n 


V„ 


V, 


n+1 


2 u 


fc —1 1 +fc fc —1 fc —1 1 +fc 


U. 


2+n' y n-l ty n 






2+n 


1+fc 1+fc 1+fc fc—1 1+fc 


l+n^2+n 17 n 


n+1^2+n 


a 'l+n a 2+n (; n <; l+n (; 2+n — u > 


or 


( 2 v k ~ 1 v k ~ 1 -I- 77 1 r 2yy fc —^^7;^ 1 -I- v 

\ za n+l u n ' Uj 2+n u 2+n) u n+ 1 \ ALL n u n- 1 ' “l+ra^l+n/ 


fc -1 


1+fc 1+fc fc+1 


2+n V 2+n u n+l 


V-n 


7y 1 + ^ ?; ^+ 1 7, 1 + fc 

^l+n^n+l tJn 


= {En _ ^-^tl+n^vlpvt 1 = Q 


i+fc„fc+i„i+fc 


‘l+n U n+l 


Wn 


We then have 


(2v k ~ 1 v k ~ l 4- u 1+k v 1+k ) i) k ~ 1 

\^ a n u n-l w u l+n u l+n) u n 

v 1+k v k+1 v 1+k 
u l+n u n+l 


= 3c, 


or 


(2n k ~ 1 v k ~ 1 J t-v k +^v k ^\ — r iriJ k ~^ 1 v kJrl yy^+l — f) 

\ zu n u n-l ' a n+l U n+l) U n oca n+l U n+l U n ~ u - 


Setting u k —> c k / 2 u k , v k —> c k ^ 2 v k , we arrive at 


( 2 ' 


.fc- 1 ,.fc -1 


— 1 I 7/* : + 1 7;k+ 1 '\ 1 _ o„.k-\-l k-\-l k-\-X — Q 

-1 i u n+l v n+l) u n oa n+l u n+l u n ~ u> 


which checks with (14.22[) . So Eas. (14.37l) - (14.40l) constitute the Lax pair for (14.2 II) - (14.221) . 


(4.43) 

(4.44) 

(4.45) 

(4.46) 

(4.47) 

(4.48) 

(4.49) 

(4.50) 


5 Conclusion 

In this paper we present a bilinear form for a integrable lattice related to the well-known relativistic Toda 
lattice. Through Hirota’s bilinear integrable discretization method a fully discrete version of the lattice is 
obtained. Bilinear BT and Lax pair for both the original and the fully discrete lattice are investigated. 
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